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Recent experiments on thin-film microcavities give evidence of Bose condensation of exciton- 
polariton states. Inspired by these observations, we consider the possibility that such exotic "half- 
light/half matter" states could be observed in thin-film organic semiconductors where the oscillator 
strength is generally stronger than in inorganic systems. Here we present a theoretical model 
and simulations of macroscopic exciton-polartiton condensates in anthracene thin films sandwiched 
within a micro-meter scale resonant cavity and establish criteria for the conditions under which 
BEC could be achieved in these systems. We consider the effect of lattice disorder on the threshold 
intensities necessary to create polartion superfluid states and conclude that even allowing for up to 
5% angular disorder of the molecules within the crystal lattice, the superfluid transition remains 
sharp. 



I. INTRODUCTION 

Over the past few years there has been consid- 
erable interest in the dynamics of cavity polaritons 
in low-dimensional inorganic microcavitiesP"S In the 
strong coupling regime, where the photon-exciton in- 
teraction is larger than the exciton and photon damp- 
ing rates, photons within the cavity and excitons 
within the material become coherently coupled to form 
bosonic quasi-particles termed "polaritons" PtED Macro- 
scopic Bose-Einstein condensate (BEC) states have re- 
cently been repo rted in GaAs and in CdTd/CdMgTe 
microcavitiea^^^' and recent reports give evidence 
of polariton condensation in one-dimensional ZnO 
microwiresP^ 

Polaritons obey a dispersion relation that splits into 
two branches with an avoided crossing near the intersec- 
tion of the cavity dispersion and the exciton resonance 
frequency. In a typical inorganic quantum-well semicon- 
ductor microcavity, the Rabi splitting between Wannier- 
Mott excitons and the cavity modes is typically on the 
order of 10 meV. However, for organic systems (with 
Frenkel excitons) exciton oscillator strengths are consid- 
erably higher than their inorganic counterparts and con- 
sequently the Rabi splittings can be on the order of 100 
me y|i5ii6] while polariton lasing has been reported in 
organic microcavities thus far, polariton BEC has not 
been in these systems. This paper is also motivated by 
the recent reports of BEC of a two-dimensional photon 
gas in a dye-filled optical microcavity^ In these exper- 
iments, the strong coupling requirement is not met and 
the quasi-particles within the cavity are indeed best de- 
scribed as photons rather than polaritons. 

In this paper, we consider the preparation and dynam- 
ics of BEC states in an organic microcavity consisting 
of an acene thin-film sandwiched within a microcavity of 



thickness L. To date, the evidence for BEC has not been 
reported in organic semiconductor microcavity systems. 
Whereas in inorganic systems, the electronic states are 
best described as Wannier excitons with weakly bound 
electron-hole pairs, in organic semiconductors the exci- 
tons are Frenkel excitons that are largely localized on 
a single molecular unit. Consequently, their electronic 
transition moments are generally larger compared to in- 
organic systems resulting stronger coupling to the photon 
field and consequently larger Rabi splittings. 




FIG. 1: Anthracene unit-cell in the {a, b} crystalographic 
plane. Arrows indicate the direction of the 5*0 — > Si elec- 
tronic transition dipole moment. 

While BEC in atomic systems involves weakly inter- 
acting particles, these systems are not uniform since the 
condensate is influenced by the form of the trapping po- 
tential. While liquid 4 He, for example, is also a BEC con- 



densate, perturbative treatments can not be used on this 
system due to strong inter-particle interactions and corre- 
lations. Consequently, polaritons offer a unique glimpse 
into BEC physics since the shape of the condensate state 
itself is unambiguous for a uniform system. The specific 
signature of a BEC state in the form of a spontaneous co- 
herence offers a unique opportunity to study localization 
effects in two-dimensional systems with varying degrees 
of spatial disorder. 

In a cavity system, the disorder can be introduced 
either by speckles within the optical field, by the ori- 
entational disorder of the molecular sites, or by site- 
energy disorder within a perfectly periodic crystal struc- 
ture. The role of disorder is an interesting issue since it 
would allow one to probe indirectly at least the transi- 
tion between the superfluid and Mott insulating regime. 
There has been some speculation that a Bose-glass phase 
is intermediate between the superfluid and Mott insulat- 
ing phases implying that disorder transforms the Mott 
insulting phase to a Bose-glass. However, other theo- 
ries suggest that disorder can change the Mott insulating 
phase directly to a superfluid. Recent experimental ev- 
idence involving trapped atoms indicate the transition 
is disorder-induced and reversible for superfluid and co- 
existing superfluid-Mott insulator phases. Furthermore, 
spatial disorder plays a key role in the fractional quantum 
Hall effect and vortex pinning in type-II superconductors. 

Our analysis begins with a molecular-level description 
of a single layer of polyacene molecules in which the local 
electronic excitations are described within the Frenkel ex- 
citon model with dipole-dipole couplings determined by 
quantum chemistry. The local excitons interact with a 
single mode of the photon field within the dipole approxi- 
mation and we determine the polariton dispersion curves 
by diagonalizing the full exciton+photon Hamiltionian 
for a given configuration of the lattice molecules. This 
treatment allows us to take an arbitrary configuration of 
the molecules within the lattice and generate polariton 
dispersions. In Sec. Ill we adopt a field-theoretical ap- 
proach developed by Carusotto and Ciuti which treats 
the cavity as a non-equilibrium steady-state systemP^SIII] 
This allows us to identify steady-state polaritons as pa- 
rameterized by the intensity of the driving field. We then 
use the Bogoliubov prescription to separate fluctuations 
from the stationary solution and use this to determine 
the excitation spectra for the lower polariton branch. 
We then consider the response of the system to inhomo- 
geneous perturbations within the cavity and show that 
above a threshold intensity for the driving field, the scat- 
tering signal collapses to a single point indicating the 
onset of superfluidity. 

The cut-off temperature for the amplification is ulti- 
mately determined by the binding energy of the exciton. 
Because this is large in organic semiconductors (hundreds 
of meV), parametric amplification should be readily ob- 
servable at significantly higher temperatures than in in- 
organic quantum-well microcavities, even at room tem- 
perature. 
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FIG. 2: Polariton/exciton dispersion curves for (a) co-parallel 
(«— • • • <—) and (b) co-facial (t ■ • ■ t) aligned dipoles. 



II. THEORETICAL MODEL 

Polaritons are composite quasi-particles that form due 
to the strong coupling between excitons and the radiation 
field. Denote a) as an operator that creates an excitation 
with energy e x and a as an operator which removes an 



excitation with {a, a}} — 1 and with similar operators 

for the radiation field, [V'fc, V^/] = ^fcfc' which create and 
remove photons with wave-vector k. The Hamiltonian 
describing the coupled exciton/photon system is given 
by 



H = e 



n, 



x qtg + tkJki'ii'k + -TrCo^fc + 4>ja) (1) 



where is the Rabi frequency which depends upon 
the number of photons in a given mode k, the oscillator 
strength of the transition, and the orientation of the exci- 
tation's transition moment relative to the polarization of 
the photon. Close to the resonance condition Huik = e x , 
the energy eigenstates of system are best described in a 
mixed or entangled representation of the system given by 



cos ffc sin tffc 
— sin 8k cos 6k 



m 



(2) 



where the mixing angle is given by tan2#fc = £ln/(e x — 
hu>k) where the U\. and L,\ operators are quasi-particle 
operators which create polaritons in the upper (U) or 
lower (L) branches with energy 




n 2 Q 2 R 



where + denotes the upper polariton (UP) and — denotes 
the lower polariton (LP) branches respectively. 

Within a microcavity, the photons obey standing-wave 
boundary conditions such that the transverse modes are 
quantized giving rise to the dispersion 



e k = (hc/? ] )^/k 2 + (2iTn/Ly 



(3) 
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FIG. 3: Polariton energy levels vs. photon wave-vector for 
a finite anthracene slab. The thick blue curve corresponds 
to the lower polariton branch (LP) while the thick red curve 
is the cavity dispersion. At two upper polariton branches 
can as well as evidenced by the avoided crossings within the 
exciton bands, (a) Exciton-polariton dispersion for perfectly 
ordered anthracene lattice. (b)Exciton-polariton dispersion 
when 88 — 15° of orientational disorder is included in the slab 
model. (c,d)Exciton f x and photon f p fraction for dispersions 
in (a,b) 



where rj is the refractive index and n = 1, 2, • • • is the 
index of the transverse mode. About k — 0, the cavity 
dispersion is quadratic 



2 1,2 



€k « A- 



2m* 



(4) 



where m* = 2irrj/cL is the effective mass of the cavity 
photon and A = 2nhc/Lrj is the cut-off energy for the 
cavity. From this point onward, operators ^ and i/»k 
create and remove photons within the microcavity and 
index k refers to the photon wavevector in the plane of 
the slab. 

We can extend this simple model to describe the inter- 
action between a photon field and slab of molecular sites 
within the cavity. First, we assume that the excitons are 
local to the molecular sites and can be described within 
the Frenkel exciton model using the HamiltoniarP 

Hex = '^e i a\a i +'^Ji j {a\aj+a i j a i ) (5) 

i i=£j 



tons at site i located at rj. In the absence of internal re- 
laxation processes, the exciton dynamics within the film 
are determined by the exchange integral Jy . Within the 
point-dipole approximation, this can be written as 



V 2 I r% 



(6) 



where fli are transition moments for the S — > Si elec- 
tronic transition on a given site and 



Tj is the 



vector connecting the centers of molecules i and j. Fig. [I] 
shows the "herringbone" arrangement of the anthracene 
molecules within a single crystallograpic plane. 

The coupling between the cavity modes and the exci- 
tons is given by 



n/2 = 7 £ 



e • Hi 



(e^aj^ + e-** 1 ^) (7) 



where N is the number of sites within the 2D slab and e 
is the polarization vector of the photon. The model used 
here has been used in a variety of contexts, notably by 
Agranovich and Gartstein in their recent study of low- 
energy exciton-polaritons in one-dimensional systemsP 
The site-model is particularly attractive since it al- 
lows us to parameterize it using both spectroscopic and 
quantum-chemical information about a specific molecu- 
lar system. 

One should, in principle, add to our model the fact that 
the excitons are composite quasi-particles consisting of 
bound electron/hole pairs. For organic systems, the ex- 
citon binding energy is in the range of 0.4 to 0.5eV and 
the electron and hole largely reside on the same molecular 
site, which justifies the use of the Frenkel exciton model 
for these systemsPsHHI An important feature which we 
do neglect, however, is the coupling between the local 
site exciton and the internal vibrational motions of the 
molecule itself. In organic systems, such vibronic cou- 
plings play important roles in the spectroscopy and dy- 
namics of conjugated systemsP^El 

Since we are interested in polaritons in molecular crys- 
tals, we need to consider how the crystallographic ar- 
rangement and possible molecular disorder will affect 
the formation of polariton states within such a system. 
With this in mind we consider an ordered 2D array of 
acene molecules within a microcavity. At 290 K, poly- 
acene systems typically form monoclinic crystals with 
two molecules per unit cell in the plane containing the 
molecular y axis oriented alternately ±20° from the crys- 
tallographic b direction. When viewed from above the 
plane, we see the "herringbone" packing geometry char- 
acteristic of polyacene molecular crystals. The electronic 
transition moment for the S a — > S± transition lies along 
the molecular y axes and hence lies in a plane parallel to 
the reflectors. 

Consider the case of a molecular dimer (such as two 
polyacenes) interacting with a photon field where we let 

t 



where the operators a] and hi create and remove exci- hi and h[ denote exciton operators for one monomer 
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and a 2 and a\ be exciton operators for the other with 
the exchange interaction denoted by J. In the "local" 
{ipk, 01,02} basis, we can write 



+ J(a\a 2 + &2&1) 



+ n fl /2(^(ai + a 2 ) + h.c) 



(8) 
(9) 



Taking e± = e 2 , H can be brought into block-diagonal 
form by creating Davydov exciton states via 



D± = (Si ± a 2 )/v 2 



(10) 



which results in 



H = ?iu k 1>l# k + (ei + J)DlD + 



(ii) 



Taking the transition moments parallel in the same plane 
where 9 is the angle between the transition moment and 
the centers of the two molecules, i.e. /*• ■ ■ ■ the ex- 
change coupling is given by 



J 



(2cos 2 (6») 



\e)) 



(12) 



For < 9 < cos^l/v 7 ^), J < and the dimers form 
a J-aggregate while for angles cos" 

_1 (1/V3) < < it/2, 
J > and the dimers form an H aggregate. At 9 — 
cos _1 (l/\/3) = 54.7°, the exchange coupling vanishes 
and no exciton splitting occurs. 

We now note that only the symmetric Davydov state 
is coupled to the photon field and we can use a unitary 
transformation to define the upper and lower polariton 
branches 



D + 

i>k 



cos <p k sin (p k 
-sin0 fc cos^fe 



(13) 



where the mixing angle is given by tan20^ = Qr/(ci — 
hujk) where the U k and L, operators are quasi-particle 
operators which create polaritons in the upper (U) or 
lower (L) branches with dispersion 



± (ei + J)+/iw fc 

h k = o T 



'(ei + J) - huj k 



^|14) 



Consequently, for a molecular dimer, we have three 
states: the upper and lower polarition states and a dark 
exciton state. The relative ordering of the exciton energy 
levels depends upon the sign of J, which in turn depends 
upon angle between the co-facial stacking planes. For J- 
aggregates, J < the bright state is lower in energy than 
the dark state where as for H-aggregates J > and the 
dark exciton state is lower in energy. In Fig. [2] we com- 
pare the dispersion curves for a molecular dimer taken as 



a J-aggregate with 9 = and an H-aggregate (9 = it/2). 
Here we use parameters: hu!k = k 2 /2m* + A, m* = 0.01, 
Q,r = 0.5| J\, and A = for purposes of this discussion. 
For the J aggregate, the lower Davydov state carries a 
net coupling to the photon mode and the lower polariton 
branch is formed solely by coupling to this state. The 
upper Davydov state does not contribute the formation 
of polaritons. The situation is exactly reversed for the 
H-aggregates in that only the upper Davydov state is 
coupled to the photon field. 

This basic idea can be expanded to describe the full 
2D slab in two ways. First, we can write the transition 
moment for the unit cell by taking the vector sum and 
differences jl± = {fL\ ± fl 2 )/V2 of the molecular transi- 
tion moments within the cell. This results in transition 
moments fi± which are orthogonal and directed along the 
a and b crystallographic axes. Thus, taking the dipole- 
dipole coupling to be short ranged so that only neigh- 
boring cells are coupled, the exciton dispersion can be 
written as 

E x (k) = e Q + 2J a cos(k x x) + 2Jb cos(k y y) 

where J a and J& are computed using J a = [i 2 _/a 3 ri and 
Jb = where a and b are the lattice constants. For 

anthracene: a = 8.5262A and b = 6.038A are the lattice 
constants for the planar array shown in Fig. 1. Since 
the wavevector k refers to the wavevector in the lattice 
plane, this dispersion is flat compared to the dispersion 
for the cavity photons and one can immediately extend 
Eq.[l4]for the 2D lattice. 

A more realistic approach is to include orientational 
and energetic disorder directly into the model and obtain 
polariton states by numerical diagonalisation. In the cal- 
culations presented next, each molecular site is described 
by a transition moment fj,i = 1. 1901a. u. corresponding to 
an oscillator strength of f osc = 0.1555 as determined by 
ab initio (CIS/6-31G*) calculations on an isolated an- 
thracene molecule using Gaussian '03)221 Rather than us- 
ing the ab initio vertical excitation energy for anthracene, 
we set the site excitation energy as ej = 2.0eV. In the 
end, this value simply establishes a reference point and 
all excitations and energies can be re-scaled relative to 
the vertical exciton energy for purposed of comparison 
with experiments. Each dipole is located at the crys- 
tallographic site and oriented alternately ±20° from the 
crystallographic b direction. This coupling gives rise to 
exciton couplings of — 96.2meV between the two acenes 
within the unit cell. We also assume a static dielectric 
constant of 77 = 3, a value typical for organic solids. Fi- 
nally, we model the cavity by taking a photon energy 
cut-off of A = 1.6eV which is detuned from the exciton 
resonance energy. Using this approach, we can easily in- 
clude the effects of energetic and orientational disorder on 
the exciton dispersion. Thus, armed with the molecular 
crystal geometry, the excitation energies, and the elec- 
tronic transition moments, one can determine the polari- 
ton energies varying with photon wave- vector for a given 
microcavity, by diagonalizing the above Hamiltonian for 
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a given photon wave-vector. 

In Fig. [3}i we show the single-particle eigenvalues of 
the coupled exciton/photon system for photon A:- vectors 
between ±25/im as determined by diagonalizing 

H = H ex + + H ex / ph 

at specific values of k. Close to k = 0, the lower polariton 
branch is nearly parabolic and one see clear evidence of 
anti-crossings. The lower thick-line highlights the "lower 
polariton" (LP) branch in this system. Looking closely, 
one can identify multiple upper polariton bands that pass 
through the density of states. Also it is important to note 
that not every exciton state couples to the photon field. 
These dark states result from the fact that their net tran- 
sition dipole is either very small or nearly perpendicular 
to the polarization vector of the driving laser field. 

Within an organic molecular crystal, orientational dis- 
order of the individual molecules inhibits the formation 
of delocalized excitonic states within the lattice. Conse- 
quently, this may limit the ability of the system to form 
collective coherent excitations such as polaritons. We 
can easily examine the effect of orientational disorder by 
randomly orienting the dipoles about some mean angle 
6 = 20° with variance a 1 . In Fig|3jD we show the po- 
lariton curves for a slab with a$ = 15° of orientational 
disorder. The primary effect of the orientational disor- 
der appears to be that polaritons formed by coupling 



to the lowest Davydov states splits into multiple polari- 
ton states as evidenced by the increased number of anti- 
crossings that can be spotted within the density of states. 
Moreover, the first anti-crossing occurs at slightly lower 
values of k. In Fig. [3](c,d) we show the relative mixing 
between the cavity photon and the anthracene slab for 
both the ordered slab (c) and the orientationally disor- 
dered slab (d). In both cases the mixing is quite strong 
even as one moves from the anti-crossing region in spite 
of the fact that the cavity cut-off is w 0.2eV off resonance 
from the exciton band. 



III. FIELD THEORY TREATMENT 

While a fully "molecular" treatment is possible in prin- 
ciple using a site-wise representation, the ~ 1/xm length- 
scales suggested by the dispersion curves in Fig. [3] would 
require two dimensional slabs with upwards of 10 6 an- 
thracene sites. On the other hand, the bare exciton dis- 
persions are essentially flat over the range of wave vec- 
tors that are important for polariton formation within 
the cavity. Bearing this in mind, let us define field op- 
erators ipx and ipe for the exciton and cavity modes and 
write the Hamiltonian density within the Hartree-Bose 
approximation: 



H = dx 



+ / dx( J B(x,t)^t(x)+ J B*(x,t)Vc(x)) (15) 



where w c and uj x generate the dispersion curves for the cavity and exciton modes and 7 C and ~/ x describe the cavity 
and exciton decay rates and Q is the Rabi operator which couples the cavity field (at a given photon k-vector) to the 
exciton field as given by Eq. [7] 

Steady state solutions for the fields can determined by first writing E(x,t) = E p (x) exp(— iu) p t) and VvW 
, : exp(— iojpt) and writing 

Lb x -ij x /2 + g\<j) x . v \ 2 /2~Lj p n/2 \ { <f> x , s (x)\ ( \ ( , 

d/2 u c llc /2 -u>J'\ <Mx) ) + { E p (x) ) a ^ ib ^ 

This leads to a set of equations for the steady-state, 

Lu x -i lx /2 + g\(f> XiV \ 2 /2-Lo p ft/2 Y\f 



•■>x,s \ _ 
t>c,s 



(1/2 Qj c -i 1c /2-u p ) \E p {x) 



(17) 



These can be solved iteratively for a given driving field, E p (x). Furthermore, if one assume that the system is 
homogeneous and the laser field consists of a single mode at k p , the steady state equations simplify to 

t> x .s\ = _fu} x 0s. p )-ij x /2 + g\(j> XiV \ 2 /2-w p ft(k p )/2 \ _1 / \ - . 

h, s ) \ fi(kp)/2 uj c (k p ) - i lc /2 - uj p J \E P J- 

where the uj x and uj c are the exciton and cavity dispersions evaluated at k„ of the driving field. 

I 



A. Lower Polariton Branch neglect any interband interactions. We thus introduce 

Since we are primarily interested in what happens close 
to k = 0, we consider only the lower branch polariton and 



field operators ^l(x) and ^(x) and write the Hamilto- 
nian (with H = 1) in the mixed representation as as 



u 



LP 



dx*i,(x)(wz, P - i7/2)^(x) 



dx(vf[(x)) 2 (vI/ L (x)) 5 



+ / dx(^(x,t)*i,(x) +ft.c.) 



(19) 



where the term i)ip(- iV) describes the lower polariton 
curve obtained from slab calculations described above. In 
doing so, we can incorporate site-wise disorder directly 
into the continuum/field theory equations. 

The field operators obey commutation relations: 
[*»(x), *J (x')] = 5(x - x')Sij. The physical system itself 
is in a driven non-equilbrium state and we need to in- 
clude a decay rate 7 that encompasses both cavity leak- 
age and non-radiative effects, and an external driving 
field (i.e. the laser pulse) E(x,t) as measured within the 
cavity. Lastly, g describes the polariton-polariton scat- 
tering. In the calculations that follow, we take these to 
be phenomenological parameters since both 7 and E(x., t) 
depend upon the particular experimental situation and 
microcavity cell. The non-linearity parameter can be 
approximated by g — Anha/m* where a is the S- wave 
scattering length. It is difficult to obtain a good esti- 
mate for this parameter so we set g = 0.1 meV which 
would correspond to a polariton/polariton collision fre- 
quency of g/h = 0.152ps _1 which seems to a reason- 
able estimate at low exciton densities. We are currently 
working on a more robust estimate of g based upon exci- 
ton/exciton scattering in molecular crystals.^ Likewise, 
we take 7 = 0.1 meV since this produces a decay time 
which is short compared to the exciton radiative lifetime 
but longer than the Rabi period of (~ 10 fs) that can be 
estimated from the exciton oscillator strength. 

Within the mean-field approximation, the time evolu- 
tion of the fields (?/H x ) = (^l(x))) is given b> |34 | 3S | 



5V. m j 

Mm. 



(20) 



This results in the following non-linear, inhomogeneous 
Schrodinger equation: 



.fly>(x) 



{Co LP - 47/2 • 



5 |V(x)| 2 Wx) 



(21) 



For the case of atomic condensates (and no dissipation 
or driving fields) this is the Gross-Pitaevskii equation. 

We again assume the material to homogeneous and 
that the laser fie ld can be described as a plane wave 
within the cavityP^ 21 * 23 ! In this case, the steady state 
equations take the form of a Ginzbcrg-Landau equation 
of state for the driven system. 

(w LP (k p ) -lu p - i 1 /2 + g\il> s \ 2 )i) s +E P = (22) 
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FIG. 4: Mean field energy vs. field strength for k p 
0.55/im- 1 . 



In general, the stationary solution ip s is complex- valued. 
In Fig. [6^, we plot the mean-field interaction energy fflV'sl 2 
vs. external field strength E p for k p = 0.55/xm _1 which is 
close to the inflection point on the LP dispersion curve. 
Here we see that below a critical value, E P}Crit , of the field 
intensity, gIV'sl 2 increases steadily as the field strength 
increases, jumps sharply before increasing linearly. 

One can consider the pumping intensity to be propor- 
tional to the number of polaritons in e Bose gas created 
by the laser field. Consequently, the sudden change in 
the mean field energy indicates the transition from a or- 
dinary Bose gas to the superfluid or condensate state. 
As we discuss in the Section [ill Bl this transition is an in- 
dication that the steady state solution is parametrically 
unstable over a continuous range of k values either side 
of k p . 



B. Excitation spectra 



Having determined the steady-state solutions of the 
homogeneous gas (Eq. 22 1, we need to determine whether 



or not such solutions are stable with regards to impos- 
ing a inhomogeneous perturbation. For this, we derive 
the normal modes of the field to fluctuations about the 
stationary solution by writing the field in terms of sta- 
tionary term and a fluctuation ip = tp s + J-0) 34 l 35 l Substi- 
tuting this into the energy functional above and writing 
rl = "Ho + Hi + H.2 + • ■ ■ . 
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H = / dxip*(x)({j LP - u p - i 7 /2)V> s (x) 



+ | / dx|^(x)| 4 + y dx(£(x)^ s *(x) + / l .c.) 

Hi = y"dx^ t (x)[(^p-^-z 7 /2 + g |^(x)| 2 )^ s (x) + ^(x)]^(x) 
+ /i.e. 

H 2 = j dx5$(*)(<b LP ~ta p - i-f/2 + 2.g|V^(x)| 2 )^(x) 

+ | f rfx((v,;(x)) 2 w(x)) 2 + (^(x)) 2 (^t (x)) ^ 



(23) 



(24) 



(25) 



Of these, Hi = for the steady-state solution and where the operator C is given by 

Ho gives the energy of the stationary state. Since H2 

is quadratic in the fluctuation operators, stability of 

the stationary state can be determined by analyzing 

the normal-mode modes given the eigenvalues of the 

Boglioubov-de Gennes equation.^ 



(£ - u p ) 



lij(x) 




U;(x) 


_ «i(x) _ 


= 6 





(26) 



£ = 



- 3 (^) 2 e- 2 ^- x 



(27) 



Since we have assumed the system to be translationally invariant, wave- vector k remains a "good" quantum number 
and we can Fourier transform the eigenvalue problem and write 



,± 



^(k) 



where 



u LP {\i)+g\i> s \ 2 -i~{/2 9i> 2 s 

-g{r a ? -w iP (2k p -k)- 5 |V> s | 2 -*7/2 



(28) 



(29) 



At each value of k, we find two eigenvalues, (^(k)*, 
which in general are complex valued. Those with posi- 
tive imaginary components are unstable modes meaning 
that non-equilibrium population within these modes will 
grow exponentially to form condensate states while those 
with negative imaginary components are stable modes 
and population will decayP^The "idler branch", w~(k), 
related by symmetry to the "signal branch" (w + (k)) via 



~(2k p — k). The excitations are are 



due to the coupling between a generic mode at k and an 
idler mode at 2k p — k induced by the driving field. This 
corresponds to the process where by two photons from 
the driving field are converted into a signal-idler pair: 
{k p , k p } — > {k, 2k p — k}. This coupling between modes of 



different k is characteristic of Bogoliubov quasi-particle 
excitations (bogolons) and parametric amplification. 

In Fig. [5] we show the eigenvalues of Ck for the condi- 
tions giving rise to Fig. [4] Here we see that the idler and 
signal branches intersect and become degenerate about 
the points of intersection. This is known as "branch 
sticking". Increasing E p causes the branch sticking to 
spread over the entire intersection region. We also see 
that 3(w + (k)) > indicating that these modes have be- 
come parametrically unstable. Increasing E p > E p ^ cr i t 
causes all modes within the intersection region to be- 
come unstable. In Fig. [6j [7] we show the range of k that 
gives rise to unstable modes. For weak field intensity, all 
modes are stable for all k. As E p increases, first one re- 
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FIG. 5: Real and imaginary eigenvalues of Ck for below (a,b) 
and above threshold (c,d) E VtCr n. (parameters: E p — 0.02eV 
(a,b), E p = 0.05eF(c,d), uj p = 1.5968eV, g = O.lmeV, k p = 
Q.b'o^m" 1 along the a;-direction,and 7 = O.lmeV) 



gion then a second region becomes unstable. These even- 
tually merge such that above the critical threshold, all 
u> + (k) within the intersection region have positive imag- 
inary eigenvalues. It is also important to notice that the 
energy dispersion is linear (or nearly so) over the entire 
branch-sticking range. This signals the on-set of a super- 
fluid state according to Landau's criterion. As a conse- 
quence, non-equilibrium population within these modes 
will grow exponentially to form condensate states. 



C. Effects of orientational disorder 

We now consider the effects of orientational disorder 
in the film on the critical behavior of the field strength 
E Pt crit and, consequently, on the stability of the system. 
The shape and position of the lower polariton curve, and 
thus E pcrit , depends on the orientation of the molecu- 
lar units inside the sample. Even at criogenic temper- 
atures, one expects some degree of disorder in the ori- 
entation of the molecular sites within the lattice. We 
can introduce disorder by sampling the orientation of the 
anthracene molecular unit about its crystallographic ori- 
entation (see^ for details) . Since both the exciton dis- 
persion and coupling to the photon field are sensitive to 
dipoles' orientation, each sampling will generate a unique 
LP dispersion. We therefore use this to compute E p ^ cr i t 



for different realizations of the lattice for a given angular 
variance erg. 

On Fig. [9|we show the distribution of E P:Cr i t for 3 dif- 

, I — I / F7771 . . „ ,,. . . . 



ferent values of erg (see^. For a perfect crystalline lattice 
erg = and E v ^ crit ~ 0.04ey. As one expects, for small 
erg E p<cr it is distributed about the E p crit and increasing 
erg results in broader distribution of E P}Cr i t . 



Fig. 10 shows 



values of E p ^ cr i t for various choices of erg, averaged over 
500 configurations. Systematic variation of the E p ^ cr i t 
with erg is small for low disorder; however, for larger erg 
the behavior of the curve changes and the slope increases. 
For a system held under cryogenic condition one expects 
erg to be very small and consequently we expect that the 
transition from a normal Bose gas to condensate should 
be sharp even with some degree of molecular disorder. 
However, increasing the lattice temperature, the transi- 
tion will be difficult to observe. 



D. Response to inhomogeneous perturbation 

In the physical system, point defects within the thin 
film or defects in the reflecting surface introduce potential 
scattering sites for the polarion field. In order to examine 
the effect of such sites we introduce a perturbing field, 
/ and consider its effect on the stationary states found 
above. Taking ijj{t) = ip s + St/j(t) as the time-dependent 
field in the presence of the scattering sites, this evolves 
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FIG. 6: Mode stability versus applied field strength. Shaded 
region indicates modes with positive imaginary frequency 
components. Above E PtCr it all modes within the branch stick- 
ing region are unstable. 
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in time according to 
. d 



i-H =C-5^ + f 



(30) 



As above, we seek stationary solutions to this equation. 
Since the steady-state part satisfies the homogeneous 
part of this equation, we can consider the (stationary) 
response to the perturbation by solving the inhomoge- 
neous part 



<ty>(x) = -C- 1 ■ /(x) 
Thus, taking /(x) = V(x)(ip s , — i/j* ), we have 



<5V(k) 



k) 



zr 1 



(31) 



(32) 



where Vk is the Fourier transform of the perturbing field. 
In essence, the polaritons are elastically scattered over 
all fc-vectors about the constant energy shell given by 
^Lp(kp). The resulting signal is then given by ^(k)! 2 . 

In practice, we create a two-dimensional array (100 x 
100) of coordinate grid-points, select randomly 300 sites 
and at those sites set V(x) = O.lmeV. At all other sites 
y(x) = 0. The inhomogeneous terms V(k) and V(k — 
2k p ) are then found by fast-Fourier transform (FFT). 
In the latter case, we use the shift theorem followed by 
FFT. In this way, we can rapidly average over multiple 
realizations of random scattering events and compute an 
averaged signal intensity. 

We consider the response for a cavity that is off- 
resonance with the Davydov exciton as depicted in the 
polariton dispersions in Figsj3] Here the cavity cutoff 
is about O.leV below the lowest Davydov exciton en- 
ergy. In Fig. [8] we show the response signal (plotted as 
log(|(5^(k)| 2 )) for the conditions given in Fig. p[c,d), i.e. 
close to the LP inflection point. The orange/yellow ring 
indicates that scattering occurs to all fc-vectors on the en- 
ergy shell given by uji J p(k s ). A polariton droplet created 




FIG. 7: Dependence of positive values of Im(a; + (k)) on the 
applied filed strength and the range of in-plane wave vectors k 
where system looses its stability. Flat region represents modes 
with negative imaginary components where steady state is 
stable (fluctuations' population decays). However there is a 
laser regime in which Im(oj + ) becomes positive and popula- 
tion with this mode grows. Here ct> p =1.5968 eV, fc p =0.55 /im, 
and E p , crit =0.039 eV. 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

FIG. 8: Scattered intensity distribution for k p — 0.55^tm _1 
(plotted as log(|<5i/>(k)| 2 )) in the Rayleigh scattering regime 
for E p below the critical intensity E PiCr u. 



under these conditions would evolve on the "on-shell" 
ring as observed in RefPl 

We next consider what happens when the driving field 
is close to the bottom of the LP dispersion curve and 
nearly resonant with the LP polariton. In Fig. 11 the 
driving field is detuned lmeV to the blue from the ujlp- 
Here, the idler and signal dispersion touch at k p creating 
a single point of degeneracy that becomes unstable as E p 
increases. 



In Fig. [12p , and b we show the signal resulting from 
the polaritonic defects. Below threshold, all eigenvalues 
of Ck have negative imaginary components and we see 
a tight ring indicating that scattering occurs to all fc- 
vectors on the Wip(fcp) constant energy shell (Rayleigh 
scattering). Above threshold intensity, this ring collapses 
to a single peak at k = {fc p , 0} indicating the formation of 
a single component "condensate" or polariton superfluid 
state. This does not occur at fc = 0, rather exactly at 
fcp. This is consistent with what one would expect for 
momentum distribution for an interacting Bose gaa^ 

?i(k) = z 5(k p ) + n(k) 

where z D = N a /N is the condensate fraction and n(k) 
is the non-condensate occupation, which is more or less 
constant with a weak divergence ~ k^ 1 as k — > kp.^l 
This state is special since it is the only one that grows as 
the number of Bosons in the system N increases. Note 
that the peak is shifted since the system is driven by the 
pumping field that has a wave- vector k p . Since we have 
assumed the pumping field to be of the form i?pe 4kp ' x , it 
equivalent to imposing a moving reference frame on the 
system. Furthermore, our theory does not include inelas- 
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FIG. 9: Probability density function of E p , cr i t for different degrees 
of disorder 69; cj p =1.597eV, k p =0.55 fJ.eV; circles, squares, and 
diamonds correspond to 50=2°, 59=5° and 69=10° respectivily. 
Each set of data is fitted with Gaussian distribution shown by solid 
line. Best fit is obtained by Gaussians with: /i = 0.0388 and a = 
0.0053 for 0-0 = 2°;^ = 0.0381 and a = 0.0130 for 0-0 = 5°; 
fi = 0.0333 and a = 0.0245 for ag = 10°. 



tic scattering events that would scattering the signal off 
the energy-shell. 



IV. DISCUSSION 

In this paper, we discuss the conditions under which 
polartiton BEC could be achieved in an organic micro- 
cavity with J-aggregate molecules such as anthracene. 
We have identified conditions under which spontaneous 
coherence and parametric amplification may be achieved 
in a polariton gas giving rise to a quantum phase transi- 
tion to a superfluid phase in an organic thin-film system. 
This opens the door to experimental probes of polari- 
ton condensate dynamics and construction of a quantum 
phase diagram for organic micro-cavity systems. To date, 
evidence for BEC in organic microcavities has not been 
presented in the literature although efforts to achieve this 
are underway. 

It is important to discuss factors not included in our 
model which may inhibit or perhaps preclude the for- 
mation of stable condensate states in these systems. Of 
primary concern is the ability to fabricate well-ordered 
molecular arrays within the microcavity. It is absolutely 
critical that the exciton transition moment be parallel to 
the reflecting walls of the cavity in order to maximize the 
coupling to the photon field. It is also likely that local 
orientational disorder will inhibit the formation of long- 
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FIG. 11: Real and imaginary eigenvalues of Ck- (parameters: 
E p = 0.02eV (a,b), lo p = u LP {k p ) + l.OmeV, g = O.lmeV, 
k p = 0.15^m _1 ,and 7 = O.lmeV) 



ranged coherences within the cavity. As seen in Fig. [3] 
orientational disorder gives rise to an increase in anti- 
crossings and an increase in the number of intermediate 
polariton bands. These may serve to deplete the LP pop- 
ulation. Likewise, strong coupling between the electronic 
and intramolecular vibrational modes leads to a detun- 
ing of the spontaneous emission band (fluorescence) from 
the vertical absorption absorption band. However, since 
the strongest electron photon coupling in these systems 
is local to a given molecular site, this will decrease as 
the Davydov state contributing to the LP is delocalized 
over multiple molecular sites. Our current approach docs 
not include internal relaxation due to inelastic scattering 
events. From detailed balance considerations, inelastic 
scattering will shift the scattering signal towards lower 
wave vectors and hence the condensate peak to shift to- 
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FIG. 12: Signal intensity in the super-fluid regime (close to 
k = 0) for E p below (a) and above (b) the critical threshold. 



wards k = 0. Thermalization must be rapid compared 
to the polariton lifetime in order to achieve a true BEC 
state at k = 0. 
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